Creation and characterization of vortex clusters in atomic Bose-Einstein condensates 
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We show that a moving obstacle, in the form of an elongated paddle, can create vortices that are 
dispersed, or induce clusters of like-signed vortices in 2D Bose-Einstein condensates. We propose 
new statistical measures of clustering based on Ripley's K-function which are suitable to the small 
size and small number of vortices in atomic condensates, which lack the huge number of length 
scales excited in larger classical and quantum turbulent fluid systems. The evolution and decay of 
clustering is analyzed using these measures. Experimentally it should prove possible to create such 
an obstacle by a laser beam and a moving optical mask. The theoretical techniques we present 
are accessible to experimentalists and extend the current methods available to induce 2D quantum 
turbulence in Bose-Einstein condensates. 
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I. INTRODUCTION 

Turbulent superfluid helium has been shovifn to exhibit 
features typical of ordinary, classical turbulence, on scales 
larger than the average inter-vortex spacing, such as the 
same Kolmogorov energy spectrum in three dimensional 
(3D) turbulence [THS] as well as differences, such as non- 
Gaussian velocity statistics [71 IMTTj . Two dimensional 
(2D) turbulence is very different from 3D turbulence. In 
3D turbulence, on scales larger than the average inter- 
vortex spacing, there is a flow of energy from large scales 
to small scales through the Richardson cascade, whereby 
large eddies are broken up into smaller and smaller eddies 
giving a Kolmogorov energy spectrum with k~^/^ scal- 
ing. In 2D turbulence there is an inverse energy cascade, 
where energy flows from the scale of energy injection 
(small scales), to larger scales as like-charged vortices 
cluster [12]. This phenomenon is established in classical 
fluids and is thought to be the mechanism responsible for 
Jupiter's great red spot [131 E]. The concept of an in- 
verse energy cascade was introduced by Onsager [15J for a 
point-vortex gas consisting of many many vortices, where 
it was found clusters of like-signed point vortices have a 
negative temperature (further detail can be found in sub- 
sequent analyses [TniHTj). Beyond the seminal work of 
Onsager, there is no fundamental theory that exists for 
small systems of vortices. In Bose-Einstein condensates, 
while there have been a number of theoretical investiga- 
tions [rSlET] and some experimental work [55] , the ques- 
tion of whether an inverse cascade is a feature of a system 
of turbulent vortices in quasi-two dimensional systems re- 
mains open; this is firstly because these condensates are 
relatively new and secondly (and more importantly) be- 
cause, due to their relatively small size, they lack the 
large range of length scales typical of other 2D turbulent 
flows (such as planetary atmospheres). 

On the other hand, ultra-cold Bose condensed atomic 
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gases are ideal to investigate the dynamics of few- vortex 
systems as well as turbulent dynamics of many vortices. 
The dimensionality of such condensates can be easily 
controlled, allowing direct study of vortices in two and 
three dimensions. Advanced experimental methods for 
the imaging and detection of vortices in Bose-Einstein 
condensates (BECs) have been developed [23] and new 
theoretical proposals have shown that vortices can be 
easily manipulated in BECs |21]- There are many tech- 
niques that can be applied to nucleate vortices in BECs, 
such as: engineering the condensate phase profile [251155] : 
stirring the condensate with a blue or red detuned laser 
(for experiments see [13 HI] and theory [5il [5^] PU]): 
mixing and merging condensates of well defined phases 
[?TH55] : moving a condensate past a defect [55]; rotat- 
ing the trapping potential or thermal cloud |34ri39j ; and 
cooling the condensate with a rapid quench through the 
phase transition (Kibble-Zurek mechanism) l40l-[42|. Fi- 
nally, vortices can be nucleated from dynamical instabili- 
ties, such as through the decay of the snake instability of 
a soliton [331 03] , the bending wave instability of a vortex 
ring [45l [46] or surface mode excitations of the conden- 
sate [37] HH] . In the first experimental demonstration of 
3D quantum turbulence in a BEC, tangled vortices were 
created by shaking the condensate with an oscillatory 
trapping potential [47H50] . 

In this paper we add to the numerous existing tech- 
niques applied to induce vortices in BECs, and apply a 
moving object with an elliptical paddle shape to create 
vortices. We demonstrate that the trajectory of the op- 
tical paddle through the condensate can be controlled to 
create vortices that are both well distributed or clustered 
into groups with like-winding. To determine if vortices 
are indeed clustered, and if this clustering increases or 
decreases significantly with time, we develop tools which 
are more suitable to the small size and the relative small 
number of vortices which can be generated in ultracold 
atom BECs. Drawing on the wealth of available statis- 
tical pattern recognition methods, we analyze our data 
using Besag's function [51], a modification of Ripley's 
K-function [55] |S3] , which has been extensively applied 
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across a variety of scientific fields to measure clustering 
and clumping of discrete objects e.g. [54H60] . Motivated 
by Besag's function, we develop some new measures of 
independent clustering when the system is comprised of 
two unique types of discrete objects. In our case, these 
are vortices with '+' or ' — ' winding in a BEC. These 
techniques can distinguish between the cases of mixed 
clusters of vortices and independent clusters of vortices 
of like- winding in the condensate. After reviewing our 
theoretical model, we analyze vortex structures obtained 
by different forms of stirring by applying both Besag's 
function and nearest neighbor techniques. Finally, we 
summarize our main conclusions, that the resulting vor- 
tex clustering is strongly dependent on the trajectory of 
the moving paddle. We find no compelling evidence of 
an inverse cascade in these small systems, in the sense 
that we don't see an increase in clustering over time. 



II. NUMERICAL MODEL 

We simulate the dynamics of a trapped two- 
dimensional Bose-Einstein condensate stirred by an opti- 
cal paddle by integrating the following 2D dimensionless 
time-dependent Gross-Pitaevskii equation (GPE): 



= _lv2v^+'^V^ + Fp7A + K2<j|V|V, (1) 



dtp 



where the interaction strength, = 2v27raiV/a2, is 
written in terms of the scattering length, a, and the 
total number of condensate atoms, N. V is scaled as 
ip = ttr"^ I^/N, where 4" is the condensate wavefunction 
and is normalized to unity i.e. / dx|-!/)p = 1. Lengths 
and time are scaled as f/ar = r and tw^ = where r 
and t are dimensional with units of meters and seconds 
respectively. Oz = ^Juj^^rmi^ and = yjh/ {mwr) are 
the axial and radial harmonic oscillator lengths, deter- 
mined by the axial and radial trapping frequencies, lo^ 
and bJr respectively. 

The potential, Vp, describes a far-off-resonance blue- 
detuned laser beam shaped into a paddle that follows a 
rotating and circular stirring trajectory given by 



Vp — Vo exp 



2 ixcosiut) — ysvLiiujt)) 


{y cos(wt) + X sm{ujtyf' 
rf2 



(2) 



where x = x ~ vsm{t) and y — y ~ wcos(t). Vq gives 
the peak strength of the potential and is selected to be 
Vb ~ 2.6/1, where fi is the chemical potential of the BEC. 
•q determines the paddle elongation and d the width. Ex- 
perimentally, the paddle can be shaped by shining a far- 
off-resonance blue detuned laser through a mask as in 
[61]. For a paddle rotating with a frequency, w, at the 
center of the condensate, x = x and y = y. For a pad- 
dle moving at a constant radius from the center of the 



condensate without rotating, we take uj = \. In all sim- 
ulations the paddle is initially linearly ramped up to its 
maximum stirring frequency, w, after which the conden- 
sate is stirred at constant u until ts = 20. The paddle 
is then ramped off linearly over t = 5 by making the 
replacement 



Vp 



(3) 



m Eq. For paddle sizes and stirring frequencies which 
result in vortex formation, the condensate dynamics are 
evolved for a further time of f = 45. 

Numerically, Eq. ([T]) is solved pseudo-spectrally with 
periodic boundary conditions and integrated in time by 
applying an adaptive 4th-5th order Runge-Kutta method 
with the help of xmds (55] . The initial state for our simu- 
lations is obtained by a short propagation of the 2D GPE 
in imaginary time (by making the replacement t ^ —it in 
Eq. ([l])), and applying a stationary paddle potential. The 
condensate is parameterized by the nonlinear interaction 
strength K2d = 10399. For example, by choosing exper- 
imentally relevant axial and radial trapping frequencies 
of uJz = 2tt X 50 Hz and a;^ = Stt x 5 Hz, this describes a 
condensate of 2.2 x 10^ ^^Na (6 x 10^ ®'^Rb) atoms with 
scattering length a = 2.75nm (a — 5.29nm). Our simula- 
tions are run on a grid of spatial extent —20 to 20 with 
gridsize Ng — 512 (see appendix). 



III. VORTEX GENERATION 

This paper looks at three stirring motions of the pad- 
dle, Vi, Vii and V^n, each of which are specific realizations 
of Eq. ([2]), corresponding to a rotating paddle, a paddle 
moving on a deferent from the condensate center and a 
paddle moving along a deferent while rotating, respec- 
tively (see figure 1, top row). The stirring motions of 
the paddle studied generate contrasting vortex configu- 
rations, as discussed below. For each case, we track the 
total number of vortices nucleated, as shown in figure 
3a. In all cases, at later times, vortices are lost to the 
edge of the condensate. The condensate edge is selected 
by identifying where the density falls to less than 30% 
of the maximum condensate density at that time. The 
resulting profile is then smoothed to give the condensate 
edge. When two vortices are closer than a critical sepa- 
ration distance and are of opposite winding, vortex pair 
annihilation occurs, a mechanism which also reduces the 
total number of vortices in the condensate. 

For all simulations the total angular momentum of the 
condensate is also tracked, given by 
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FIG. 1: (Color online) Column 1: Density slices for paddle Vi with parameters: d = 0.5, oj = 6 (dimensionless units); Columns 
2-4: Density and vortex position images for paddles Vi, Vii and Viii. Parameters: Vi (column 2): d = 1, lj = 4; Vu (column 
3): d = 1, ljj = 1, V — 4; Viu (column 4): d — 1, uj = 2, v = 4. All paddles have Vo = 150, 77 = 8 (dimensionless units). Positive 
vortices are identified by (pink) + signs and negative vortices by (blue) o signs. Time is indicated in white at the bottom left 
hand side of each plot. 



Case I: Paddle rotating at the condensate center 

Firstly we look at a paddle rotating about its center 
with frequency w at the center of the condensate. This 
is modeled by evolving Eq. ([I]) with Vi — Vp{x ^ x, y ^ 
y). The smallest paddle size we consider is d = 0.5, 
rotating at a frequency w = 6. The rotating motion of 
the paddle produces circular spiral sound waves, which at 
late times interfere with each other giving a wave interfer- 
ence pattern (see the first column of figure 1). Paddles 
with a larger width {d = 1), rotating at frequencies of 
cj = 4, 6 and 8, nucleate vortices in addition to creating 
spiral sound waves. The density profile for a paddle ro- 
tating at frequency a; = 4 is show in the second column 
of figure 1 (see supplementary movie 1 for a comparison 
of paddles rotating at frequencies w = 4 and w = 8). A 



greater rate of rotation increases the number of vortices 
initially nucleated, as expected. In all cases, vortices are 
initially nucleated from the ends of the paddle with wind- 
ing opposite to the direction of rotation of the paddle, as 
depicted in figure 2. At subsequent times, when the local 
superfluid velocity surpasses the critical velocity for vor- 
tex nucleation |63j, vortices of both signs are nucleated 
from both the center and ends of the paddle. A cen- 
tered rotating paddle imparts a small amount of angular 
momentum to the condensate. The angular momentum 
imparted is proportional to the frequency of the rotat- 
ing paddle (see figure 3b). Note that after paddles of 
frequency a; = 6 and 8 have been ramped off, the con- 
densate angular momentum saturates to approximately 
the same value. 
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FIG. 2: (Color online) Schematic diagram of vortices nucle- 
ated from a paddle rotating at the condensate center, Vi as 
the paddle is ramped up to w = 8. Positive vortices are identi- 
fied by (pink) + signs and negative vortices by (blue) o signs, 
the paddle profile is shown in purple. 



Case II: Paddle stirring at constant radius 

The second stirring motion of the laser paddle we sim- 
ulate is a paddle stirring the condensate at a constant 
radius from the condensate center, or traveling along a 
deferent, Vu = Vp{u} = 1) (see figure 1, column 3, top 
row, for a schematic of the paddle motion). This paddle 
trajectory creates vortices which tend to group initially 
in like-signed clusters as shown by the progressive time 
samples of the condensate density profiles in figure 1, col- 
umn 3 (see also supplementary movie 2). The conden- 
sate gains a large amount of angular momentum (refer 
to figure 3b), and consequently at final times there is a 
significant imbalance in the total number of vortices with 
positive and negative winding. 




(b)Angular Momentum 

FIG. 3: (Color online) Variation of a) total vortex number 
and d) condensate angular momentum with time. Vortices 
induced by a rotating paddle, Vi, with d = 1 and v = 4,6 and 
8 are represented by orange pluses, purple crosses and pink 
diamonds respectively. Data corresponding to the paddles 
with V = 4 and d = 1, moving with trajectories Vu (cj — 1) 
and Vui {lj — 2), are given by green triangles and blue squares 
respectively. 



Case III: Paddle rotating and stirring at constant radius 

The final stirring trajectory of the laser paddle we sim- 
ulate, Vui = Vp, is a combination of the two previous 
motions, with the paddle moved at a constant radius 
from the condensate center while rotating at a small fre- 
quency. This motion can also be described as a paddle 
rotating and moving along the deferent (see figure 1, col- 
umn 4, top row, for a schematic of the paddle motion). 
The effect of adding the rotational motion of the paddle 
to its trajectory mixes the clusters produced, resulting 
in smaller groups consisting of 3 and 4 vortices (com- 
pare the vortex distributions in figure 1 columns 3 and 4 
and supplementary movie 2). While initially (for times 
t < 12) the angular momentum transfered to the conden- 
sate is reduced in comparison to that of paddle Vu (see 
figure 3b), this stirring trajectory results in the greatest 
transfer of angular momentum to the condensate at later 
times. 



IV. CLUSTERING ANALYSIS 

A. Ripley's K-function 

We analyze the clustering of vortices formed by stirring 
a 2D condensate with a paddle by applying Ripley's K- 
function, a statistical pattern analysis method used as a 
measure of spatial clustering. In the context of clustering 
of like-signed vortices, Ripley's K-function is dependent 
on the total number of liked signed vortices, N , within 
the total condensate area. A, and can be expressed as 

. N N 

i=i j=i 

where fij{r) ~ 1 for a vortex, j, within a distance r 
of the reference vortex, z, with like-winding. Otherwise, 
fij{r) = if i = J, or if the distance between vortex i 
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FIG. 4: The function fij{r) in Ripley's K-function. 



and j is greater than r. That is 



1 Vry <r,i^j 
V Vij > r or i ~ j 



(6) 



here rij is the distance from a reference vortex i to the 
comparison vortex j with hke-winding. This is depicted 
in figure 4. 
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FIG. 5: Besag's function L{r /tc) (see Eq. 9) for vortices of 
negative winding at varying times (see legend) for a conden- 
sate with vortices nucleated by a paddle with trajectory Vu. 
Parameters: d = 1, a; = 1 and v = A. 

Ripley's K-function looks at the number of like-signed 
vortices within a radius, r, from the position of an arbi- 
trarily chosen vortex, i, at its center (see figure 4 and 5). 
If the number of vortices with like-winding per unit area 
within this radius, r, is greater than the overall number 
of like-signed vortices per unit area for the whole con- 
densate, then the vortices are said to be clustered. Clus- 
tering results in K{r) increasing faster than if vortices of 
either sign are distributed in a spatially random manner, 
that is, if they follow a Poisson distribution. Ripley's-K 
function for a poisson-distributed data set takes the form 



K{r) = 



For a linear scaling of Poisson-distributed 



data, it is useful to normalize Ripley's K-function to 
H{r) — yjK{r)/-K. Ripley's L-function, also commonly 
known as Besag's function, is obtained from further nor- 
malization of Ripley's K-function: 



L{r) = ^K{r)/TT-r. 



(7) 
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FIG. 6: (Color Online) Schematic diagrams of independent 
(A) and co-clustering (B) in systems with two distinct types 
of objects, represented by pink pluses and blue circles respec- 
tively. 



As the condensate area. A, does not necessarily remain 
constant over all times, we scale r by the characteristic 
condensate radius Vc = \J Ajix for that time and in our 
subsequent analysis evaluate 



L(r/r,) 



A 



N N 



which simplifies to: 



i=i j=i 



r 



(8) 



L{r/rc 



N N 



iV2 



i=l i = l 



r 



(9) 



Besag's function is zero for like-signed vortices which are 
randomly distributed, takes positive values for vortices 
clustered over that spatial scale, and is negative if the 
vortex distribution is dispersed. That is: 



L{r/r,) = 



1 Clustered 
Random 
-1 Dispersed 



(10) 



The radius around a centered vortex containing, on aver- 
age, the most like-signed vortices per area, is called the 
radius of maximal aggregation, and is given by the value 
of r which maximizes L{r) |59) . 

For a paddle rotating at a constant radius from the 
condensate center, evaluating L{r/rc) (see equation 9) 
for positive vortices, as seen in figure 5, shows that the 
clustering of vortices decreases with time. Although the 
vortices are clustered, the amount of clustering is not 
constant or increasing in time, giving no evidence of an 
inverse cascade for this system. 



B. Measures of independent clustering of 
like-signed vortices 

While Besag's function gives a measure of the cluster- 
ing of vortices with the same winding, it does not dis- 
criminate between cases where like and opposite signed 
vortices are clustered in the same spatial region and cases 
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FIG. 7: (Color online) Comparison of evolution second near- 
est neighbor data, C2{t), (top figure) and fourth nearest 
neighbor data, C4{t), (bottom figure) for paddle trajectories 
Vn, Viii and Vi with v — 8 represented by green triangles, 
blue squares and pink diamonds respectively. Time is mea- 
sured from the beginning of stirring the condensate. 



FIG. 8: (Color online) Comparison of evolution second near- 
est neighbor data, C2{t) (top figure), and fourth nearest 
neighbor data, C4{t) (bottom figure), for vortices created by 
rotating a paddle at the condensate center {Vi) with v = 4 
(orange +), v — 6 (purple crosses) and i; = 8 (pink dia- 
monds). Time is measured from the beginning of stirring the 
condensate. 



where like-signed vortices are clustered in spatially inde- 
pendent regions (refer to figure 6 for a schematic illus- 
tration). It is necessary to make this distinction when 
looking for a measure of an inverse cascade process, as 
the clustering of like signed vortices is not expected to 
occur in the same spatial region as clustered vortices 
of the opposite sign. To address this issue we define a 
new measure of clustering which uses the sign of nearest 
neighboring vortices to determine if clustering occurs in 
regions that are spatially independent. 

We express this measure of clustering based on looking 
at the sign of all j neighboring vortices up to the i?th 
nearest neighbor of an arbitrary reference vortex i as 

^ N B , . 

Here c^- = 1 if the vortex i and its j"' nearest neigh- 
bor are of the same sign and Cij = if vortex i and its 
J*'' nearest neighbor are of opposite sign. If vortex i is 
separated by distance greater than Rc = rc/3 to its j"* 
nearest neighbor, then Cij =0. _B is the maximum near- 
est neighbor to the reference vortex i. It is necessary that 



the value of Rc chosen is greater than the average inter- 
vortex separation distance and on order of the largest 
cluster size. Vortices closer than R^. to the condensate 
edge will bias the calculation of Csit) as they have an 
area less than ttR'^ surrounding them, which their nearest 
neighboring vortices could inhabit. To correct for these 
edge effects we omit these vortices which are less than 
Rc from the condensate edge from the set of reference 
vortices, but still include them in the set of comparison 
vortices for vortices a distance greater than Rc from the 
condensate edge. Systems for which Cb — 0.5 are ran- 
domly distributed and when C b takes values greater than 
0.5 the objects are clustered. We note that these mea- 
sures can be applied generically to investigate cases of co 
and independent clustering of two discrete objects and 
could be simply extended to look at co and independent 
clustering of many discrete objects. 

A comparison of the evolution of CB{t) for the simula- 
tion runs described previously is shown in figures 7 and 
8. From the nearest neighbor analysis we learn: 

• From figure 7, we can see that a paddle moving at 
a constant radius from the condensate center (Vii) 
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creates vortices that are initially very clustered, in- 
dicated by C2 and C4 taking values very near one. 
After the paddle is turned off the vortices remain 
clustered, with C2 and C4 not decreasing below 0.5 
until t « 50. At long times {t > 50) vortices be- 
come randomly distributed. 

• A paddle moving with trajectory Vm creates clus- 
ters that are initially smaller in size than purely 
moving the paddle at constant radius from the con- 
densate center {Vu). After the paddle is turned 
off, the vortex distribution closely follows that of 
run Vu with vortices remaining clustered until long 
times (t > 50) when they become randomly dis- 
tributed. 

• When the paddle is only rotated at the condensate 
center, Vi, vortices never become clustered (refer to 
figure 8). 

For all cases, regardless of how the vortices are initially 
nucleated, evaluating Csit) gives no evidence of a ten- 
dency of increasing clustering of like-signed vortices over 
a scale of Rc in a turbulent 2D BEC after the laser paddle 
has been ramped off. Our observations imply that clus- 
ters of like-signed vortices exist due to the way in which 
they were induced in the condensate but do not naturally 
tend to cluster. 



V. CONCLUSION 

In this paper we have covered two main objectives: 

• We have extended the available methods for cre- 
ating vortices in 2D atomic Bose-Einstein conden- 
sates, demonstrating that a paddle can be used to 
stir a condensate in two quite different ways, cre- 
ating long-lived vortex clusters or more randomly 
distributed vortices that are turbulent in two di- 
mensions. 

• A new statistical measure of clustering based on an- 
alyzing nearest neighbor vortices was defined, moti- 
vated by a well known statistical spatial point pat- 
tern analysis technique, Besag's function. These 
measures have been applied to analyze how vortices 
are distributed in 2D condensates. 

We find that a paddle moved through the condensate 
at a constant radius from the center creates vortices of 
both positive and negative winding in clusters. When the 
paddle is rotated at the condensate center, vortices cre- 
ated are initially clustered co-dependently in the same 
local spatial regions and later disperse throughout the 
condensate. For a combination of both moving a pad- 
dle at a constant radius through the condensate while 
simultaneously rotating the paddle, the vortices induced 
are less clustered than if the paddle is only moved at a 
constant radius from the condensate center. The later 



method can be applied to create long-lived vortex clus- 
ters in BECs. The clusters are considered long lived in 
terms of two relevant sets of timescales; the constraining 
timescale determined by realistic experimental lifetimes 
of BECs and the timescales determined by the system 
size and properties. The relevant timescales intrinsic to 
the system size and length-scales are determined by the 
average separation distance between vortices (Isep) and 
is given by Tsep = ^lepll where 7 is the phase winding of 
a vortex and the largest turnover time t = 7rr^/27r « 35, 
is determined by the condensate size. As the largest 
turnover time is smaller than the longevity of clusters, 
which persist for ^ 50 (see C(4) in Figure 7) we describe 
the clusters to be 'long-lived'. For our choice of experi- 
mental parameters, cluster lifetime is ~ Is, which is also 
long-lived in comparison to the typical experimental life- 
time of BECs, (- 10s). 

The extent of clustering was quantitatively measured 
by evaluating two statistical measures of clustering; ap- 
plying a modified Ripley's function and a technique based 
on comparing the sign of nearest neighboring vortices. 
We did not observe an increase in clustering over time. 
This was despite evolution times longer than the largest 
turnover time r = 7rr^/27r « 35, determined by the con- 
densate size. Our system contains too few vortices to 
determine if the relevant physical process for 2D turbu- 
lent systems in atomic Bose-gases is an inverse cascade of 
incompressible kinetic energy from small to large scales 
manifesting in a clustering of like-signed vortices. In par- 
ticular, it would be difficult to apply traditional methods 
used for large systems (planetary atmospheres, super- 
fiuid helium) , based on Fourier transforming the velocity 
field and analyzing the spectra of energy and enstrophy 
over many decades in wavenumber space. Our statisti- 
cal analysis based on Ripley's K-function and on nearest 
neighbor methods provides a way to quantify an increase 
or decrease in the degree of vortex clustering. As these 
methods are constructed from a knowledge of the position 
and winding of vortices in the system they are readily ac- 
cessible experimentally. Information on vortex location 
in condensates can be obtained experimentally through 
standard absorption imaging techniques, e.g. [23, 34, 35], 
and winding of vortices is found by analysis techniques 
giving phase information, such as condensate interferom- 
etry, e.g. HZHMUSS]. 
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Appendix A: Error Checking 
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(a)Total Vortex Number 



To check our results are independent of gridsize used, 
the simulation grid size was doubled from Ng = 512 to 
Ng = 1024 for a paddle rotating at the condensate center 
with V — 6 and d — 1. In figure 9 a comparison is made of 
the total vortex number and condensate angular momen- 
tum, when the grid size is doubled. The angular momen- 
tum is calculated by evaluating Eq Q. The reasonable 
agreement between rates of vortex production and elim- 
ination, as well as evolution of the condensate angular 
momentum in both runs establishes that the gridsize of 
Ng = 512 applied in the simulations presented in the 
body of the paper is adequate. The small variance in 
results from doubling the gridsize are attributed to the 
condensate edge selection routine used. A further source 
of difference is the chaotic nature of vortex dynamics in 
turbulent systems. A small amount of numerical noise 
would be enough to seed a difference in vortex trajecto- 
ries. 



O- N, = 512 




(b)Angular Momentum 



FIG. 9: b) Evolution of angular momentum and a) total num- 
ber of vortices nucleated by a paddle rotating at the center of 
the condensate (Vi), with v = 6, d = 1. Simulation gridsizes 
Ng = 512 (purple) and Ng = 1024 (pink). 
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